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We present a theoretical study of plasmonic lattice solitons (PLSs) formed in 
two-dimensional (2D) arrays of metallic nanowires embedded into a nonlinear 
medium with Kerr nonlinearity. We analyze two classes of 2D PLSs families, 
namely, fundamental and vortical PLSs in both focusing and defocusing 
media. Their existence, stability, and subwavelength spatial confinement are 
studied in detail. © 2010 Optical Society of America 

OCIS codes: 190.6135, 240.6680. 

When the size of photonic systems is reduced to subwavelength scale, the light confinement 
and propagation are severely limited by diffraction. Overcoming this major challenge of con- 
fining and manipulating optical energy at the nanoscale has attracted a rapidly growing 
interest in nanophotonics [1]. To this end, one promising approach employs surface plasmon 
polaritons (SPPs) of metallo-dielectric nanostructures [2]. Although most of the experimen- 
tal and theoretical studies of SPPs have focused on their linear optical features, there is 
a growing interest in the optical properties of SPPs in the nonlinear regime [3-7]. This is 
primarily because the strong enhancement of the optical field induced by the excitation of 
SPPs can significantly boost nonlinear optical eflFects, as well as the promising applications of 
nonlinear SPPs in the all-optical control at nanoscale. A new approach to achieve subwave- 
length confinement and control of the optical field, employing one-dimensional (ID) arrays 
of parallel metallic nanowires embedded in an optical medium with Kerr nonlinearity, was 



1 



recently proposed [6]. Under certain conditions the optical nonlinearity induced by SPPs 
field exactly balances the discrete diffraction in the plasmonic array, and as such collective 
excitation of ID plasmonic lattice solitons (PLSs), are formed in the array. Because the ra- 
dius a of the nanowires and their separation distance rf, are much smaller than the operating 
wavelength A, the spatial width of the PLSs can be significantly smaller than A. 

In this Letter, we study the nonlinear optical modes in 2D arrays of metallic nanowires 
which are embedded into a Kerr medium [Fig. 1(a)]. We found that 2D PLSs of subwave- 
length spatial extent in the transverse section could exist in this geometry. In addition to 
the fundamental PLSs, 2D plasmonic arrays support vortical solitons, i.e., solitons with an 
internal topological phase change of 2ti along a closed contour around soliton's phase sin- 
gularity. The vortical PLSs can also be of subwavelength extent. Fundamental and vortical 
PLSs are found to exist in both focusing and defocusing media. 

Our analysis is an extension of the coupled mode theory developed in Ref. [6] for ID 
coupled metallic nanowires. Briefly, the total electric and magnetic fields are expanded as a 
superposition of the modes of a single nanowire (assumed to have only the fundamental TM 
mode), and using the conjugated form of the Lorentz reciprocity theorem one obtains 

i ^^'^ + ^(0m,n+l + 0m,n-l + 0m+l,n + 0m-l,n) + 
/i((/)^+l^^+l + (/) 
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where z is the longitudinal coordinate, (})m^n is the normalized mode amplitude in nanowire 
(m^n)^ K and ji are the coupling coefficients between neighboring and next-neighboring 
nanowires, respectively, and 7 is the normalized nonlinearity coefficient. Note that 7 > 
(7 < 0) for a medium with self- focusing (self-defocusing) Kerr response. We emphasize that, 
unlike in the case of conventional dielectric waveguides, for our plasmonic array, ^ < 0, 
/X < 0, so that the linear dispersion, kz = 2tv[cos{kxd) + cos{kyd)] + Aiicos{kxd)cos{kyd) ^ 
forms a concave surface with anomalous diffraction occurring at the F point [k^ = ky = 0)^ 
while normal diffraction at the M points {kx = ky = ivr/rf) [Fig. 1(b)]. 

The soliton solution to Eq. (1) are sought in the form (j)m,n{z) = '^m,n exp(i/3/^z), where 
the amplitudes Uj^^n are independent of z and /3 is the soliton wave number (normalized by 
hi). The soliton is characterized by its power P = X^^^ |0m,nP- This ansatz is inserted into 
Eq. (1) and the resulting system of equations is solved numerically by a standard relaxation 
method. The corresponding amplitudes are then used to reconstruct the fields of the PLSs. To 
quantify the transverse size of solitons we use the effective radius, R = {J |Ep[(x — Xo)^ + 
{y-yof]dxdy/ JJ^^\E\^dxdyy/^, where (xo,?/o) = JJ^^{x,y)\E\^dxdy/ JJ^^\E\^dxdy is 
the mean center position. In the analysis below, we assumed that the nanowires are made 
of Ag and used the Drude model, 6^(cl;) = 1 ( l- ) ^ permittivity of the metal. For 
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Ag, ujp = 13.7 X 10-^^ rad/s and u = 2.7 x 10-^^ rad/s. For simplicity, we assume that metal 
is lossless. 

We first search for fundamental solitons. As in the case of the ID PLSs [6], we also find 
two types of 2D PLSs, i.e., staggered and unstaggered solitons. In the case of unstaggered 
(staggered) PLSs the phase difference of the mode amplitude in adjacent nanowires is equal to 
zero (tt). The spatial profile of the amplitude and longitudinal field of unstaggered (staggered) 
PLSs corresponding to a nonlinear change of 6n^\ = —0.05 {Sn^\ = 0.05) are shown in 
Fig. 2(a) and (b) [Fig. 2(c) and (d)]. Importantly, due to the inverted linear dispersion 
relation, staggered (unstaggered) solitons are formed in self-focusing (self-defocusing) media, 
which is opposite to the case of dielectric waveguide arrays. Note also that the solitons have 
subwavelength size, the radius being R = 0.35A [Fig. 2(a,b)] and R = 0.50A [Fig. 2(c,d)]. The 
dependence of soliton radius on /3 is shown in Fig. 2(e). The plot shows that both staggered 
and unstaggered PLSs shrink monotonically as their wave number moves away from the band 
edges. This is because staggered (unstaggered) PLSs originate from the (extended) Bloch 
modes at the point M (F), and thus when their wave number is close to the band edge, 
solitons have a large width. A larger change of nonlinear refractive index tunes the soliton 
mode deeper into the gap [see inset of Fig. 2(e)] and thus PLSs are increasingly smaller 
than A. Importantly, our calculations show that the power carried by such PLSs cannot be 
smaller than a threshold value, as can be seen from Fig. 2(f). The nonmonotonic relationship 
P = P{f3) has an important implication on soliton stability, as will be shown below. 

In addition to the fundamental PLSs, 2D plasmonic arrays also support a new type of 
PLSs, which has no counterpart in ID systems, i.e., vortex PLSs; these belong to the class 
of discrete vortices [8]. Here we focus on off-site vortices, i.e., compact vortex states whose 
singularity is located between lattice sites. Similarly to fundamental PLSs, vortex PLSs exist 
both in self-focusing [Fig. 3(a,b)] and self-defocusing media [Fig. 3(c,d)]. As can be easily 
seen from the phase profiles, the vortex has a topological charge equal to 1 (2D PLSs with 
topological charge 2 also exist, although not shown in here). Interestingly, the phase profile 
for PLSs in focusing media features a "staggered" pattern at the sites far away from the 
soliton center, which is a typical feature of "gap" vortices [9]. However, one should mention 
that the staggered vortex PLS resides at the semi-infinite gap (above the band) and thus 
is not a "gap" vortex. This is again a consequence of the inverted linear dispersions. The 
properties of vortex PLSs are summarized in Fig. 3(e,f). Several differences from fundamental 
PLSs are observed. First, unlike the fundamental PLSs where their value of /3 can be very 
close to the band edges, vortex soliton can only exist in the gaps at a finite distance from the 
band edges. This is due to the fact that the vortex mode is not a Bloch mode of the linear 
plasmonic array. Secondly, as vortex PLSs possess four intensity maxima, vortex solitons 
always have a larger width as compared to that of fundamental ones (under the same change 
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of nonlinear refractive index). Therefore, forming a subwavelength vortex generally requires 
stronger nonlinearity. Nevertheless, we find that subwavelength vortex PLSs can be achieved 
under experimentally accessible conditions. For example, the vortex presented in Fig. 3(a,b) 
has radius R = 0.3A requiring a nonlinear change of refractive index Sn^\ = 0.11. Note 
that an index change of ^0.14 was reported in Ref. [10]. Further, the PLSs size can be 
significantly reduced if the operating wavelength is scaled down, thus the requirement for 
a strong nonlinear change of the refractive index can be relaxed. Finally, we mention that 
as in the case of fundamental solitons, vortex PLSs also feature power thresholds for their 
formation. 

A relevant issue associated with the 2D PLSs is their stability. To address this issue, 
we integrate Eq. (1) by using the Runge-Kutta method with the input condition being 
the stationary soliton solution plus random noise. Our extensive simulations show that the 
fundamental PLSs obey the Vakhitov-Kolokolov stability criterion namely, staggered (un- 
staggered) PLSs are stable when dP/d/3 > {dP/d/3 < 0) and are unstable otherwise. In the 
stability region, soliton maintains its original shape over a sufficiently long distance while 
the initial added noise is rapidly eliminated. In contrast, unstable solitons tend to relax 
to stable ones with the same power and large oscillations are observed during propagation 
[Fig. 4(a,b,e)]. The vortex can be also stable, although its stability region is more limited. 
An unstable propagation of vortex PLS is shown in Fig. 4(c,d). The vortex quickly loses its 
screwed phase structure and decays into a fundamental PLS. Note that the stability domain 
of vortex PLSs increases at a shorter operating wavelength. 

In conclusion, we presented the first study of 2D plamonic lattice solitons in arrays of 
metallic nanowires embedded into a host medium with Kerr nonlinearity. Stable subwave- 
lengh fundamental and vortex PLSs are found to exist in both focusing and defocusing media. 
The work of N. C. P. was supported by the EPSRC. 
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/3 = —4.89, 5n^\ = —0.05. Radius R (e) and power P (f) vs /3. The solid and 
dashed lines stand for d = 7a and d = 8a, respectively. In inset, 5n^\ vs /3 for 

d = 8a, in focusing medium 8 

3 (Color online) The normalized electric field amplitude (a,c) and longitudi- 
nal component (b,d) in focusing (top panels) and defocusing (middle panels) 
nonlinear media. (a,b) /3 = 3.82, 6n^\ = 0.11. (c,d) /3 = —5.2, 6n^\ = —0.18. 
Radius R (e) and power P (f) vs normalized eigenvalue /3. The solid and 
dashed lines stand for d = 7a and d = 8a, respectively. In inset, Sn^\ vs /3 for 
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4 (Color online) Unstable propagation of the fundamental soliton at /3 = 3.8 
(top panels) and vortex at /3 = 3.82 (middle panels), showing the intensity 
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Fig. 1. (Color online) (a) Schematics of a 2D square array of nanowires with radius a and 
separation distance d. (b) The band predicted by the coupled mode theory. 
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Fig. 2. (Color online) The normalized electric field amplitude (a,c) and longitudinal com- 
ponent (b,d) in focusing (top panels) and defocusing (middle panels) media. A = 1550 nm, 
a = 40 nm, = 8a. (a,b) f3 = 3.58, Sn^i = 0.05. (c,d) /3 = -4.89, Sn^i = -0.05. Radius R (e) 
and power P (f) vs /3. The solid and dashed lines stand for d = 7a and d = 8a^ respectively. 
In inset, 5n^\ vs /3 for = 8a, in focusing medium. 
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Fig. 3. (Color online) The normalized electric field amplitude (a,c) and longitudinal compo- 
nent (b,d) in focusing (top panels) and defocusing (middle panels) nonlinear media. (a,b) 
/3 = 3.82, 5n^\ = 0.11. (c,d) /3 = —5.2, 5n^\ = —0.18. Radius R (e) and power P (f) vs nor- 
malized eigenvalue /3. The solid and dashed lines stand for d = 7a and d = 8a, respectively. 
In inset, 5n^\ vs /3 for = 8a, in focusing medium. 
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Fig. 4. (Color online) Unstable propagation of the fundamental soliton at /3 = 3.8 (top 
panels) and vortex at /3 = 3.82 (middle panels), showing the intensity profile at z = (a,c), 
z = 60 /xm (b), and z = 90 /xm (d). (e) Variation of maximum amplitude vs propagation 
distance for fundamental (dashed) and vortex (solid) solitons. 
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